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Zoom-whirl behavior has the reputation of being a rare phenomenon. The concern has been 
that gravitational radiation would drain angular momentum so rapidly that generic orbits would 
circularize before zoom- whirl behavior could play out, and only rare highly tuned orbits would retain 
their imprint. Using full numerical relativity, we catch zoom- whirl behavior despite dissipation. The 
larger the mass ratio, the longer the pair can spend in orbit before merging and therefore the more 
zooms and whirls seen. Larger spins also enhance zoom-whirliness. An important implication is that 
these eccentric orbits can merge during a whirl phase, before enough angular momentum has been 
lost to truly circularize the orbit. Waveforms will be modulated by the harmonics of zoom- whirls, 
showing quiet phases during zooms and louder glitches during whirls. 



Kepler's laws describe closed elliptical planetary mo- 
tions. A small relativistic correction accounts for the 
tiny, anomalous precession of Mercury's perihelion. If 
the sun is replaced by a black hole, the geodesic mo- 
tions can zoom and whirl in an extreme form of preces- 
sion - whirling around the center of mass in nearly cir- 
cular inspiral before zooming out along elliptical leaves. 
Zoom- whirl behavior is characteristic of strong relativity 
and could potentially be detected in the harmonics of the 
gravitational waves generated. 

There are astrophysical settings that couldpopulate 
eccentric merges, such as dense galactic nuclei [l[ or glob- 
ular clusters [2], y] ■ Consequently, it is an important as- 
trophysical question to ask: Can zoom-whirl behavior, 
an intrinsically eccentric phenomenon, survive the dissi- 
pative drain of gravitational radiation? In this article we 
report on results of numerical relativity that show zoom- 
whirl behavior in comparable mass binaries, answering 
this question in the affirmative. 

As there is no analytic description of the curved space- 
time around two black holes, we rely either on analytic 
approximations or on numerical relativity to describe 
comparable black hole pairs. Zoom-whirl behavior has 
been studied in extreme-mass ratio inspirals [3, [H, HI, @] 
and was recently found in an analytic approximation, 
specifically conservative Post-Newtonian (PN) approxi- 
mations to black hole binaries [1, Q . Now we find zoom- 
whirl orbits in full numerical relativity of spinning pairs. 
Zoom- whirl behavior has already been found in numerical 
relativity for equal-mass, nonspinning binaries in [To| . In 
that work, the initial conditions were carefully tuned to 
find a special orbit, the separatrix between bound orbital 
motion and plunge. The separatrix is studied in detail in 
the PN approximation in [9( and an analytic solution for 
the separatrix in Kerr systems was found in [III [l2| . 

In this Letter we show that zoom-whirl behavior in 
spinning pairs is a common feature of eccentric orbits 
BSE!, despite dissipation. In particular, zoom- whirl 
orbits happen well away from the separatrices and so do 
not in general require fine tuning of initial conditions 



(see also [HE3). 

Due to the computational expense of 
running these simulations, a full scan of parameters is not 
possible. To focus our investigation, we rely on analytic 
approximations to estimate good initial conditions and 
we then run full numerical simulations to easily locate 
zoom- whirl inspirals. The further utility of the analytic 
estimates is the transparency of interpretation. 

The anatomy of zoom-whirl behavior was quantified 
in Ref. [l5j where it was shown that every orbit of a 
given L can be described by one number that specifies 
the precession of the orbit per radial cycle from apastron 
to apastron. The amount by which an orbit will precess, 
that is, overshoot the previous apastron is 
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A perfectly periodic orbit looks like a closed 1-lcaf clover 
or 2-leaf clover or 3-leaf clover, or z-leaf clover. And 
each periodic orbit corresponds to a rational q made up 
of w integer number of nearly circular whirls close to 
perihelion per leaf in the z-leaf clover. The vertex v is 
more subtle and indicates the order in which the z-leaves 
arc traced out. So, a simple 3-lcaf clover for instance 
is a q — 1/3 (w = 0, v = 1, z = 3) and precesses past 
the previous apastra by A</> = 2irq = 2ir/3 per radial 
cycle. A 3-leaf clover that skips a leaf in the pattern 
each time corresponds to q — 2/3 (w = 0, v = 2, z = 3) 
and precesses by Acj) — 2nq — An/ 3 per radial cycle. 
Therefore, q quantifies zoom-whirl behavior, the integer 
part signals the whirls per leaf and the fractional part 
signals the number and order of the leaves. 

Another way to interpret the number is as the ratio 
of frequencies, q = o^/o-v, where lo^ is the average of 
the angular frequency per radial cycle and uj r = 2ir/T r 
where T r is the time between apastra. For periodic or- 
bits the frequencies are rationally related and the orbit 
will eventually close. Whirls accumulate near perihelion 
simply because the angular velocity is greatest on closest 
approach and the circumference smallest. 

Generic orbits are not periodic and do not correspond 
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FIG. 1: The 3PN effective potential. Top: L IBC o ~ 4.4. 
The straight line corresponds to the orbit r a = 30. Bottom: 
L ~ 4.1, for which the orbit with r a — 30M is the separatrix. 

to rational q. However, any generic orbit can be approx- 
imated by a nearby periodic, just as any irrational num- 
ber can be approximated by a nearby rational number. 
Kepler's ellipse corresponds to q = since it does not 
precess at all. Mercury's precessing orbit corresponds to 
q ~ 10~ 7 since it precesses very little. Technically, zoom- 
whirl behavior corresponds to q > 1, so there is at least 
one whirl, although we will generally be interested in any 
substantial precession, say q > 1/4. 

If angular momentum is fixed, q decreases monotoni- 
cally with decreasing energy. However angular momen- 
tum decreases, along with energy, as gravitational radia- 
tion is emitted and the evolution of q will not necessarily 
be monotonic. In the case of comparable mass black hole 
pairs, q will change quickly due to the rapid losses to 
gravitational waves as evident in the simulations. 

The range of zoom- whirl behavior is most easily tar- 
geted with an effective potential picture. In Ref. [3, [l5[ 
an effective potential method was used to describe the 
center-of-mass motion of black hole binaries in the con- 
servative 3PN Hamiltonian description with spin-orbit 
coupling included. There it was shown that at the turn- 
ing points of the orbital motion, the Hamiltonian itself is 
an effective potential and the motion can be read off as 
easily as interpreting a ball rolling along hills. 

Consider the effective potential on the top left of Fig . 
[JJfrom the conservative 3PN Hamiltonian [3.[l5l.[l6lll7j]. 
There is clearly a minimum of that potential which de- 
fines a stable circular orbit. An orbit energetically above 
the stable circle will execute elliptical precessions as 
shown on the top right. Note these precessions are much 
more extreme than Mercury's with q > i. It looks like 
a precession around a three leaf clover. The aperiodic 
orbit will eventually fill out an annulus. 

Also evident is a maximum of the potential, an un- 
stable circular orbit. For this angular momentum the 
unstable circular orbit is marginally, energetically bound 



since V e s just skims zero there, and has been called an 
innermost bound circular orbit (IBCO). An orbit at rest 
at infinity would asymptotically approach this circular 
orbit. And though initially of eccentricity one, e = 1, 
this orbit whirls an infinite number of times as it climbs 
the potential toward the unstable circle at the top. This 
separatrix is the infinite whirl orbit, q = oo, and is also 
known as a homoclinic orbit of eccentricity one. 

For angular momenta below this critical value, the un- 
stable circular orbit marches down in energy. For each 
unstable circle, there is a corresponding separatrix of de- 
creasing eccentricity. As an example, on the bottom of 
Fig. [TJ the separatrix between bound motion and plunge 
has apastron r a = 30M and asymptotically approaches 
the unstable circular orbit in the infinite whirl limit. So 
this is the q = oo orbit for L ~ 4.1, where through- 
out we measure angular momentum in units of \xM and 
M = mi +Ti2 while /i = 7711TO2/M . The homoclinic sep- 
aratrix is plotted alongside the corresponding effective 
potential on the bottom right of Fig. [TJ 

We expect to see zoom-whirl behavior until the unsta- 
ble circular orbit and the stable circular orbit merge at 
the ISCO (innermost stable circular orbit), the separa- 
trix with e — 0. Quasi-circular inspiral plunges at the 
ISCO and so this inflection point in the effective poten- 
tial has received preferential attention. However, orbits 
that maintain an eccentricity during inspiral will merge 
by rolling over the top of the potential, through a homo- 
clinic separatrix, behavior we observe in our simulations. 

Roughly then, for given external parameters 
(jnijm\, Si, S2), zoom-whirl behavior should be 
sought in the range Ljsco < L < Ljbco- One more 
initial condition needs to be specified to define an orbit 
and that could be either the energy or the apastron. 
We'll choose to fix the apastron. We'll fill in the details 
for a fiducial example and then flush out more general 
trends we have observed. 

The procedure is the following: (1) choose external pa- 
rameters namely mass ratio and spins (wig/mi, Si, Sj), 
(2) render the effective potential for the 3PN Hamilto- 
nian with spin-orbit couping, (3) glean the range of L 
for which there will be zoom-whirl behavior and (4) run 
simulations for this range of initial conditions. 

For our fiducial example we take mass ratio 7712/TO1 = 
1/3 and spin magnitudes S\jm\ = Sijm\ = 0.3 with 
spins antialigncd with the orbital angular momentum. 
Antialigned spins, like aligned spins, will retain the or- 
bital motion in the equatorial plane making it easier to 
distinguish whirl precession. Other than the choice to re- 
tain equatorial motion for transparency, there is nothing 
special about the fiducial configuration. We consider a 
set of orbits all of which begin at apastron r a = 30M. 
Since the simulations are so costly, we tighten the range 
to the more realistic values of angular momenta below 
Libco an d above the value of the angular momentum 
at which r — r a is the apastron of a homoclinic orbit. 
Specifically, we look between the values represented by 
the top and bottom of Fig. [JJ 4.1 < L < 4.4. 
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FIG. 2: The initial apastron is r a = 30M. Left column, 
from the top, frame 1: L = 4.10 and the orbit rapidly moves 
through zoom- whirl cycles precessing by > n (q ~ 1/2) be- 
tween apastra. The orbit rapidly loses eccentricity but still 
merges with non-zero eccentricity. Frame 2: L — 3.95 and the 
orbit precesses by nearly> 2ir (q ~ 1) in the first radial cycle 
before merging. Frame 3: L = 3.915 and the orbit is very 
nearly homoclinic, precessing by ~ 4n (q ~ 2) around the 
unstable circular orbit before radiative losses cause merger. 
Right Column: The gravitational waveforms versus simula- 
tion time were extracted at a finite radius of 60M. 

There are a few reasons why the actual range for nu- 
merical relativity will be offset from this 3PN prediction. 
For one, dissipation ensures that L changes as the orbit 
evolves. This is equivalent to the effective potential drop- 
ping as the orbit evolves. For another, the 3PN approxi- 
mation is by definition not exact and its predictions are 
not to be overstated in this strong-field regime. Finally, 
although spin-orbit coupling is included in our analysis, 
spin-spin coupling terms are not included since they in- 
troduce angular dependencies that spoil the effective po- 
tential [l5[ . Since spin-spin terms are small, their ef- 
fect will presumably result in a small perturbation to our 
3PN effective potential. Our range is intended only as a 
guide and apparently does well enough, as we'll see. 

In between this range, all orbits with our external pa- 
rameters (m2/mi = 1/3, Si/m\ = S\/m\ = 0.3), an- 
tialigned, and an apastron of r a = 30-M, will show zoom- 
whirl behavior with 1/3 < q < 00. This is the story 
told by the conservative dynamics. To test the survival 
of zoom- whirl orbits under the dissipative effects of radi- 
ation reaction we turn to full numerical relativity now. 

The runs in Fig. [2] were simulated using the Geor- 
gia Tech MayaKranc code that uses the same compu- 



tational infrastructure and methodology as in previous 
studies [H, [H, H3, HHj H3] ' namely a Baumgarte-Shapiro- 
Shibata-Nakamura code with moving puncture gauge 
conditions [23|, [24[ using the Kranc code generator (25| . 
The black-hole encounters are initiated with Bowen-York 
initial data [26[. The black holes are located on the x- 
axis: bh± are located at x + = 7.5M and X- — — 22.5M 
where m + = 3M/4 and m_ = M/4. The spatial finite 
differencing is sixth order. We used eleven levels of re- 
finement with Carpet [27) ]. a mesh refinement package 
in Cactus (28|. The finest resolution is M/143 and the 
outer boundary is at 287M. The total initial orbital an- 
gular momentum, is varied between the values of 3.9 and 
4.1. The range is shifted from the conservative PN value 
as expected, yet zoom-whirl behavior was comfortably 
found given that initial prediction. 

Fig. [5] shows three different orbits and their waveforms 
corresponding to different initial values of L. The largest 
value of L shows zoom-whirl behavior characterized by 
q ~ 1/2 before merging. The middle value of L has a 
q ~ 1 while the smallest value of L we show is very close 
to the separatrix and whirls nearly twice. We caution 
that we are only loosely reading off these values since q 
is changing so rapidly during inspiral. 

As expected, the pairs merge by rolling over the top of 
the effective potential (Fig. [I]). We emphasize that the 
nearly circular pattern of the whirl phase is not equiva- 
lent to full circularization of the orbit. In other words, 
quite importantly, the orbits shown in Fig.[2]do not merge 
through the ISCO but rather roll over the top of the po- 
tential, merging through a whirl. To compare with the 
language in a previous paper [29l |. some of those eccen- 
tric orbits were circular because they merged through a 
whirl phase and not because they merged through the 
ISCO. This could be relevant for initial conditions for 
the ringdown phase. 

The waveforms show distinctly quiet phases during the 
highly elliptical zooms followed by louder glitches dur- 
ing the nearly circular whirls. The distinctive spikes of 
zoom-whirls arc directly related to their rational num- 
ber q. This suggests that zoom-whirl orbits could be 
dug out of the data using algorithms suited for burst 
searches, perhaps in conjunction with more targeted tem- 
plate searches. 

Fig. E] plots the relative orbital separation of the black 
holes versus phase for our three cases of L. The plot 
illustrates the rapid merger of the black holes as L de- 
creases. The figure also shows that inspiral ends and 
plunge begins in the simulations for radii ~ 5M. In line 
with the numerical results, the effective potential picture 
in Fig. [T]prcdicts whirls around pcriastron, ~ 5M , which 
is much less than the PN predicted value of the ISCO, 
TiSCO ~ 8.8M, for a mass ratio of 1/3 and these spins. 
The figure therefore confirms that the zoom-whirl pair 
merge near a whirl phase and not near the ISCO. The 
top line, corresponding to L = 4.1, shows the three zooms 
before merger that we see in frame 1 of Fig. O The bot- 
tom line, corresponding to L = 3.915 and frame 3 of Fig. 
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FIG. 3: The relative orbital separation in units of M plotted 
versus the orbital phase for the series of L in Fig. [2] 

[5] is almost homoclinic with no zooms, precessing around 
the unstable circular orbit before merging. 

In general, we conclude from the PN approximation 
that zoom-whirl behavior has better odds at surviving 
dissipation and leaving a mark in the gravitational wave- 
form for (1) more disparate values of m^/mi because of 



the slower dissipation time, (2) larger spin magnitudes, 
all else being equal, because of the greater frame drag- 
ging effect, (3) larger initial apastra because of the longer 
orbital lifetime. There is less zoom-whirl behavior in 
equal-mass nonspinning pairs, hence the importance of 
sticking closer to the separatrices [T3|. Yet more gener- 
ally, zoom-whirl orbits ought to extend down below the 
separatrices. In the future, we intend to extend this work 
by pushing the numerical simulations to large mass ratio 
and extracting a measure of q. 

For any zoom-whirl pair, the glitchy waveforms fol- 
lowed by longer quiescent phases are highly distinctive 
and beg for tailored search algorithms. Furthermore, the 
loud glitches could nudge these highly processing orbits 
into an optimistic position for early direct detection. 
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